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We study a modified version of the interacting monomer-dimer (IMD) model that has infinitely 
many absorbing (IMA) states. Unlike all other previously studied models with IMA states, the 
absorbing states can be divided into two equivalent groups which are dynamically separated infinitely 
far apart. Monte Carlo simulations show that this model belongs to the directed Ising universality 
class like the ordinary IMD model with two equivalent absorbing states. This model is the first 
model with IMA states which does not belong to the directed percolation (DP) universality class. 
The DP universality class can be restored in two ways, i.e., by connecting the two equivalent groups 
dynamically or by introducing a symmetry-breaking field between the two groups. 
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A wide variety of nonequilibrium systems with a sin- 
gle trapped (absorbing) state display a continuous phase 
transition from an active phase into an absorbing phase, 
which belongs to the directed percolation (DP) universal- 
ity class |l]-|4|] . Recently, systems with multiple absorbing 
states have been investigated extensively. The interact- 
ing monomer-dimer (IMD) model introduced by one of 
us H is one of many models that have two equivalent 
absorbing states These models belong to a differ- 

ent universality class from DP. By the analogy to the 
equilibrium Ising model that has two equivalent ground 
states, this new class is called as the directed Ising (DI) 
universality class (To). When the (Ising) symmetry be- 
tween the absorbing states is broken in the sense that one 
of the absorbing states is probabilistically preferable, the 
system goes back to the DP class |n]]. Hence, the sym- 
metry between the absorbing states is the key factor in 
determining the universality class of models with several 
absorbing states. Unfortunately, no models with higher 
symmetries than the Ising symmetry (like the three-state 
Potts symmetry) are found to have a stable absorbing 
phase as yet. 

In contrast, systems with infinitely many absorbing 
(IMA) states are far less understood. All IMA systems 
studied so far belong to the DP universality class p^P^I ■ 
The number of absorbing states of these IMA systems 
grows exponentially with system size but there is no 
clear-cut symmetry among absorbing states. Recently, 
it was argued that the IMA models should belong to the 
DP class unless they possess any extra symmetry among 
absorbing states |l^,|l5|. However, no IMA model with 
an additional symmetry has been studied up to date and 
the role of the symmetry in the IMA systems is still un- 
clear. 

In this Letter, we introduce an IMA model with the 
Ising symmetry between two groups of absorbing states. 
These two groups of absorbing states are equivalent and 



dynamically separated infinitely far apart. In other 
words, an absorbing state in one group can not be reached 
from any absorbing state in the other group by a finite 
number of successive local changes. There is no infi- 
nite dynamic barrier among absorbing states inside each 
group. This dynamic barrier is similar to the free energy 
barrier between ground states of equilibrium systems that 
exhibit spontaneous symmetry breaking in the ordered 
phase. Our numerical simulations show that this model 
belongs to the DI universality class. Furthermore, we find 
that this model crosses over to the DP class by allowing 
that the two absorbing groups are connected dynami- 
cally and/or by introducing a symmetry-breaking field 
to make one absorbing group probabilistically preferable 
to the other. 

Our model is a modified version of the ordinary IMD 
model that we call the IMA-IMD model. Dynamic rules 
of the IMA-IMD model are almost the same as those of 
the IMD model with infinitely strong repulsion between 
the same species in one dimension |p[. A monomer (A) 
cannot adsorb at a nearest-neighbor site of an already- 
occupied monomer (restricted vacancy) but adsorbs at 
a free vacant site with no adjacent monomer-occupied 
sites. Similarly, a dimer (-B2) cannot adsorb at a pair 
of restricted vacancies (B in nearest-neighbor sites) but 
adsorbs at a pair of free vacancies. There are no nearest- 
neighbor restrictions in adsorbing particles of different 
species. Only the adsorption-limited reactions are con- 
sidered. Adsorbed dimers dissociate and a nearest neigh- 
bor of the adsorbed A and B particles reacts, forms the 
AB product, and desorbs the catalytic surface immedi- 
ately. Differentiation between the IMA-IMD model and 
the IMD model comes in when there is an A adsorption 
attempt at a vacant site between an adsorbed A and an 
adsorbed B. In the IMD model, we allow the A to ad- 
sorb and react with the neighboring B, so there are two 
equivalent absorbing states comprised of only monomers 
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at alternating sites, i.e., (ylOAO • • ■) and (OAOA • • •) where 
"0" represents a vacancy. In the IMA-IMD model, this 
process is disallowed. Then, any configuration can be 
an absorbing state if there are no nearest neighbor pair 
of vacancies and no single vacany between two B parti- 
cles, e.g., (• • • B0A0BB0A0A ■ ■ •). To impose the Ising 
symmetry between the absorbing states, we introduce 
the probability s of spontaneous desorption of a near- 
est neighbor pair of adsorbed B particles. At finite s, 
an absorbing configuration cannot have this BB pair. 
Hence only those configurations that have particles at 
alternating sites and no two B's at consecutive alternat- 
ing sites become absorbing states, e.g., (^40^40-BOAO • • •) 
and (0A0A0B0A ■ ■ ■). The absorbing states are divided 
into two groups with particles occupied at odd- and even- 
numbered sites (O group and E group). The number of 
absorbing states in each group grows exponentially with 
system size and there is a one-to-one mapping between 
absorbing states in two groups. It is clear that one can 
not reach from an absorbing state in one group to an ab- 
sorbing state in the other group by a finite number of suc- 
cessive local changes. Any interface (active region) be- 
tween two absorbing states in the different groups never 
disappears by itself in a finite amount of time, so there 
is an infinite dynamic barrier between the two groups. 
These interfaces annihilate pairwise only. 

The order parameter characterizing the absorbing 
phase transition is the density of active sites or kinks 
(domain walls). In the IMD model, the dimer density 
served well as the order parameter, but it cannot do in 
this model. We use the kink density as the order param- 
eter. Kinks are defined such that all absorbing configura- 
tions have no kinks but any local change of the absorbing 
configurations should produce kinks. In this model, one 
should examine, at least, three adjacent sites to check 
the existence of kinks. There are 13 possible configura- 
tions for three adjacent sites. We assign a kink to eight 
different configurations; 000, 00A, A00, BOO, 00B, BOB, 
BBO, and OBB. Five others, AOA, AOB, BOA, 0A0, and 
050, do not have a kink. In this kink representation, 
there is no mod(2) conservation of the total number of 
kinks. 

Three independent critical exponents characterize the 
critical behavior near the absorbing transition: the order 
parameter exponent (3, correlation length exponent v±, 
and relaxation time exponent i/y 0]. Elementary scaling 
theory combined with the finite size scaling theory filf l 
predicts that the kink density p(p c , L) at criticality in the 
(quasi)steady state scales with system size L as 

p(p c ,L)~L-^. (1) 

One can also expect the short time behavior of the kink 
density as p{p c , t) ~ t~ /3 / ,y ii and the characteristic time 
scales with system size as r(p c , L) ~ L v \^ v± . 

In Monte Carlo simulations, a monomer is attempted 
to adsorb at a randomly chosen site with probability 



(1 — s)p and a dimer with probability (1 — s)(l — p). 
With probability s, a randomly chosen nearest neighbor 
pair of adsorbed -B's (if there is any) is desorbed from 
the lattice. We choose the dimer desorption probabil- 
ity s = 0.5 and run stationary Monte Carlo simulations 
starting with an empty lattice with size L = 2° up to 
2 11 . The system reaches a quasisteady state first and 
stays for a reasonably long time before finally entering 
into an absorbing state. We measure the kink density 
in the quasisteady state and average over many survived 
samples. The number of samples varies from 2 x 10 5 for 
L = 2 5 to 2 x 10 3 for L = 2 11 . The number of time steps 
ranges from 10 3 to 2 x 10 5 . 

From Eq. (|l|), we expect the ratio of the critical kink 
densities for two successive system sizes p{L/2)/p(L) = 
2 /3 / !y - L , ignoring corrections to scaling. This ratio con- 
verges to unity in the active phase (p < p c ) and to 2 
in the absorbing phase (p > p c ) in the limit L — > oo. 
We plot the logarithm of this ratio divided by In 2 as a 
function of p for L = 2 6 , • • • , 2 11 in Fig. 1. The cross- 
ing points between lines for two successive sizes converge 
slowly due to strong corrections to scaling. In the limit 
L — > oo, we estimate the crossing points converge to the 
point at p c = 0.425(4) and /3/v± = 0.49(3). The value of 
(3/v± agrees well with the standard DI value 0.50. 

In Fig. 2, we show the time dependence of the criti- 
cal kink densities p(p c ,t) for various system sizes with 
p c = 0.425. From the slope of p(p c ,t) we estimate 
/3/v\\ = 0.275(5). Insets show the size dependence of the 
relaxation time r{p c ,L) and the steady-state kink den- 
sity p{p c ,L) at criticality. We estimate v\\/v±^ = 1.74(4) 
and P/v±_ = 0.494(6), respectively. All of these results 
are in excellent agreement with the DI values. 

We run dynamic Monte Carlo simulations with various 
initial configurations and get a more precise estimate of 
the critical probability p c — 0.425(1). Our estimates for 
the dynamic scaling exponents are S + r\ = 0.28(1) and 
z = 1.14(1) jL7|, where S + r\ characterizes the growth of 
the number of kinks averaged over survived samples and 
z the spreading of the active region Q . These values are 
also in excellent agreement with the DI values. 

To check the importance of the Ising symmetry among 
the absorbing states, we introduce a symmetry breaking 
field such that the monomer adsorption attempt at an 
even- numbered site is rejected with probability h [p"T|| . 
For finite h, the O group of absorbing states is proba- 
bilistically preferable to the E group. We set ft, = 0.1 and 
run stationary Monte Carlo simulations for lattice sizes 
L = 2 5 up to 2 9 . In Fig. 3, we plot ln[p(L/2)/p(L)}/ In 2 
versus p, from which we estimate p c — 0.304(2) and 
f3/v± = 0.24(1). The value of (3/v± is clearly different 
from the DI value but agrees well with the standard DP 
value 0.2524(5). More detailed study including dynamic 
Monte Carlo simulations confirms that the systems with 
finite h belong to the DP universality class p7| . 
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Similar to the case of the ordinary IMD model, 
the symmetry-breaking field makes the system behave 
like having only one (preferred) group of absorbing 
states |l(J. Evolutions of the critical interfaces (active 
region) (a) for the symmetric case (h = 0) and (6) for 
the asymmetric case (h = 0.1) are shown in Fig. 4. In 
the symmetric case, the interfaces between the O and 
E group of absorbing states diffuse until they meet and 
form a loop to disappear, which is the essential charac- 
teristic of the DI universality class. In the asymmetric 
case, the absorbing region of the unpreferred (E) group 
quickly vanishes and the interfaces between the differ- 
ent groups become irrelvant. The interfaces inside the 
preferred (O) group, which can disappear by themselves 
without forming loops, become dominant and force the 
system into the DP universality class. 

When the desorption process of a nearest neigh- 
bor BB pair is forbidden (s = 0), the system can 
find many more absorbing states with BB pairs, e.g., 
(• • • B0A0BB0A0A ■ ■ •), in addition to the two groups of 
the absorbing states for s ^ 0. These new extra ab- 
sorbing states are generically mixtures of the O and E 
group of the absorbing states. The O and E groups are 
now connected dynamically via new mixture-type ab- 
sorbing states. Consider a configuration with an inter- 
face between two absorbing states in the different groups, 
(• • ■ B0A0000A0A - ■ •), where the interface is placed in 
two central vacancies 00. With nonzero s, this configu- 
ration never evolves into an absorbing state. However, 
in the case of s = 0, it can evolve into a mixture- type 
absorbing state by adsorbing a dimer BB in the center. 
Actually, any interface can disappear by itself in a finite 
amount of time, so there is no infinite dynamic barrier 
between absorbing states. Therefore the evolution of the 
interfaces resembles the asymmetric case in Fig. 4. 

Absorbing states for s = no longer possess the clear- 
cut global symmetry which drives the system into the DI 
class. So we expect the system falls into the DP class 
like the other IMA models without an extra symmetry. 
We run dynamic Monte Carlo simulations starting with a 
lattice occupied by monomers at alternating sites except 
at the central vacant site, (■ ■ ■ A0A000A0A- ■ ■), where 
represents a defect. In Fig. 5, we plot three effective 
exponents against time; S(t), r)(t), and z(t) M. Off crit- 
icality, these plots show some curvatures. The values 
of the dynamic scaling exponents can be extracted by 
taking the asymptotic values of the effective exponents 
at criticality. From Fig. 5, we estimate p c = 0.105(1), 
S = 0.02(1), rj = 0.48(5), and z = 1.33(5). The values 
of S + r\ and z arc in good agreement with those of the 
DP values |l8| ]. Introduction of the symmetry breaking 
field h only changes the location of p c . Stationary Monte 
Carlo simulations also confirm our results Jl7[ | - 

In summary, we found the first IMA model that does 
not belong to the DP class, but belong to the DI class. 
This can be achieved by imposing a global Ising symme- 



try on the absorbing states, i.e., making the two equiv- 
alent group of IMA states that are dynamically sepa- 
rated infinitely far apart. When the symmetry between 
these groups is broken, one group of absorbing states be- 
comes completely obsolete and the evolution morphology 
changes from a loop-like to a tree-like structure, which 
ensures the system in the conventional DP class. We also 
found that the system goes back to the DP class if the 
mixture-type absorbing states between the two groups 
are added. These extra absorbing states connect the two 
separated groups dynamically and make the loop- forming 
process of the interfaces irrelvant. The absorbing states 
in all other previously studied IMA models are dynam- 
ically connected in the sense as mentioned above. This 
may explain why those models belong to the DP class. 
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FIG. 3. Plots of ln[p(L/2)/p(L)]/ln2 versus p for the 
asymmetric case with h = 0.1. 
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FIG. 1. Plots of \n[p(L/2)/p(L)]/ In 2 versus p for the sym- 
metric case. 
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FIG. 2. The time dependence of the kink density at 
p c = 0.425. The straight line is of slope 0.275 (= f3/v\\). 
Insets show the size dependence of the relaxation time r and 
the steady-state kink density p at criticality. The solid lines 
are of slope 1.74 (= V\\/v±) and -0.494 (= -p/v±). 



FIG. 4. Evolutions of the critical interfaces for (a) the 
symmetric case and (b) the asymmetric case. Monomers (^4) 
are represented by black, dimer particles (B) by grey, and 
vacancies by white pixels. 
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FIG. 5. Plots of the effective exponents against 10/t for s — and h = 0. Five curves from top to botton in each panel 
correspond to p = 0.100, 0.104, 0.105. 0.107, and 0.110. 
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